Effect of one-ion L-S coupling on magnetic properties of a correlated spin-orbit system 
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By introducing a basis for a novel realization of the SU(4) Lie algebra, we exactly solve a spin- 
orbital chain with one-ion L-S couphng (OILSC) via the Bethe ansatz (BA) approach. In the 
context of different Lande g factors of the spin and orbital sectors, the OILSC results in rich and 
novel quantum phase transitions. Some accurate analytical expressions for the critical fields are 
obtained. Both spin ordering and orbital ordering are found in a gapped singlet phase. The system 
exhibits many interesting phenomena such as nonvanishing magnetization of the singlet phase, multi- 
entrance of the singlet in the ground state when the field varies, unsymmetric magnetization in 
two-component phase and magnetization crossover for different OILSC. 



Strongly-correlated electron systems in the presence 
of orbital degeneracy have attracted much interest due 
to experimental advances in transition-metal compounds 
[1] , and the orbital ordering as well as the orbital density 
wave has been observed in a family of manganites [2]. 
Among the experimental findings are materials related 
to spin-orbital systems in one dimension, such as tetrahis 
(dimethylamino) ethylene-Cgo [3], artificial quantum dot 
arrays [4], Na2Ti2Sb20 and NaV205 degenerate chains 
[5] , and so forth. A well-known spin-orbital model is the 
SU(4) model [6] which is exactly solvable in one dimen- 
sion [7,8]. The competition of orbital degree of freedom 
with the spin results in three quantum phase transitions 
(QPT's) in a magnetic field, according to both numerical 
[9] and analytic [10] analysis. Deviations from the SU(4) 
symmetry can be caused by variations of different-site 
interaction [3,5,11,12]. Another possibility for deviation 
was considered for an SU(2)(X)SU(2) Ising-type same-site 
anisotropy in Ref. [13,10], which may occur when the 
external field turns the spin oriented in parallel to the 
orbital angular momentum frozen in some direction by a 
crystalline field. A more generally-existing same-site in- 
teraction is SU(3)®U(1) one-ion L-S coupfing (OILSC), 
the strength of which varies substantially in a range of 
0~10'*cm~^(lcm~^ '--^lK'--^10~^eV) for different elements 
[14]. As it is well understood now, L-S coupling can 
cause magnetic anisotropy [15], as far as a single ion is 
concerned. An interesting problem is to investigate the 
influence of such an OILSC on the collective properties, 
especially the QPT's of strongly-correlated spin-orbital 
systems in a magnetic field. It can be expected that the 
OILSC will bring about novel physics in the competition 
of spin and orbital degrees of freedom. 

In the present paper we shall exactly solve a spin- 
orbital chain with such an OILSC via the Bethe ansatz 
(BA) approach. In the context of the different g fac- 
tors of spin and orbit sectors, the OILSC leads to a rich 
variety of novel QPT's, with the analytic critical fields 



explicitly obtained. Both spin ordering and orbital or- 
dering exist in a gapped singlet phase, accompanied by a 
nonvanishing magnetization. The energy of the singlet is 
no longer invariant in the field due to the different g fac- 
tors, which brings about five consecutive QPT's. It also 
contributes to unsymmetric changes of the total magne- 
tization within a same phase in weak and in strong fields. 
We consider an A^-site chain with the Hamiltonian 

N 

n = Ho + nL-s + M, Wo = ^/'m+i, 

i=l 

Hl-s ^ \Y.^r ■ M^-gsH^{gssl+gtlt), (1) 

i i 

where s and I are spin-1/2 operators for spin and or- 
bit and gs and gt arc the corresponding Lande g factors, 
Pi,i+i = (2si-s,;+i + ^)(2/i-Zi+i-|-i) exhanges the neigher- 
site states. H is the magnetic field and periodic boundary 
conditions are applied throughout. Hl-s is the OILSC 
from the N sites and note that electrons have positive A 
whereas holes have negative A [15]. In the absence of the 
field, Hl-s breaks the system symmetry into SU(3)®U(1) 
when the bulk part Ho is SU(4) invariant [6]. The basis 
consists of singlet and triplet of the SU(2) Lie algebra 
{s+ -f Z+,s^ + 1^,3^^ + P}, which is valid in the pres- 
ence of the field only if g^ — gt as in the integrable spin 
ladder [16,17]. However, Hl-s and M do not commute 
due to the different g factors of spin and orbit, the con- 
ventional singlet and triplet are no longer their common 
eigenstates and fail to be the solution in the spin-orbital 
system. Nevertheless, we can diagonalize Hl-s + M a,s a 
whole, which requires a new basis for all sites 

01-y"V2 , '/'I +2/02 , 

= , „ , Vl = 03, </?2 = , „ , V3 = 04, 

where 0i = [Ti), 02 = jit), 03 = ITT), 04 = jii) de- 
note site state |s^P) and y^^ = ±g-H/X -f j{H/X), 
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^{h) = (1 + g'^h'^y^^, g± = gs ± gi- The new reahza- 
tion of the SU(4) Lie algebra takes SmnP>i — Sni^m, 

Sw=p{y)r{y)[yl+ [1/2 + s^) - ,s+(P + 1/2)], 

S30 = piy)r{y)[ys- {1/2 - F) - T (1/2 - s')], 

S12 = r{y)[l+{l/2 + s') + ys+{l/2 + P)], 

^32 = r(y)[s-(l/2 - n + yl-{l/2 - s')], S^s = s+l+, 

S02 = P{yy{y)[yh+l- - s-l+ + y{s^ - P)], (2) 



where r{y) = l/-y/l + y^, pijj) = y/ \y\, other opera- 
tors are their conjugates Snm = S^n- ^'^^ ^ certain 
choice of the basis order, the three Cartan operators Ik 

{k = 1,2,3) can be generated by commutations, e.g., 
h = [823,3,2]., h = [812,821] and h = [5'oi,^io] for 
basis order {^Q,pi,(p2,'P3)'^ ■ The new basis {po} and 
{1^1,^2,^3} form the singlet and triplet of a new SU(2) 
Lie algebra {A^ ,A~ ,A^} where A'^ = + P and 

A± = ^^2r{y)[{y - 1) {s^V - l^s^) + ^(s± + Z±)], 

though the field breaks the symmetry into four U(l)'s. 

Wc notice that a time-dependent single L-8 coupled 
particle was solved by an SU(2) gauge transformation 
[18], but here we consider a many-body system with N 
correlated ions. Note that Pij is still the permutation op- 
erator in the new basis, and Tih-S + -M. commutes with 
the SU(4) bulk part Hq. Based on the above new SU(4) 
realization we can exactly solve the model via the BA 
approach [8] with the BA equations and the eigenenergy 
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j,m 



m=l 



m— 1 



m=l 



M 



(1) 



where "^xilJ^j'^) = [l^f^ - fJ-m - xi)/{fj,f^ - + xi), 
nf^ = 0,M(0) = A/',M(4) = 0, and 1 < fc < 3; a„(/Lt) = 



^Tl^l^Vi, Eo = -J ~ |7(i7/A), El = ^ - \g+H, 



E3 = \ + \g+H and E^ = -\ + ^-/{H/X). Ni is the 
total site number in state pi. The basis order is chosen 

as {(pp-^ifip^Lpp^ifip^)'^ , where Pi E {1,2.3,4} and pp-^ is 
energetically the most favorable while (pp^ is the least 
favorable. The relation Np. = M^*"^) - M^') gives an- 
other form of energy E = J^^^i 9^^\n) + £f^^^M(2) + 
g(3)M('^). The field shifts the order from {'Po,^pi,pi2,P'3Y' 
to {p}i,ipo,ip2,p}3f for A > 0, and {ipi,ip2,^p3,^pof to 
{^1,^2,^0,^3)^ for A < 0, at point Hr = g+ |A| /{2gsgt). 
Following Ref. [19] and Rcf. [20] one can obtain the 
ground state (GS) equations for three dressed energy e^*-*, 

eW = - a2 * e'^'^- + ai * (e^^-^)" + e'^'+'l-), (3) 

where e*^°-* = e*^^-* =0 and the symbol * denotes the con- 
volution. Fermi seas filled by negative e^''~ form the GS. 



In the absence of the field, the GS is composed of all 
the four basis components when Ac- < A < Ac+ where 
Ac+ = 4, Ac- = -2ao/3 and ao = (v^7r-31n3)/2. But a 
strong positive OILSC A > Ac+ keeps the triplet gapful, 
the GS only consists of singlet and the system exhibits 
spin and orbital ordering, as shown later on. The first 
QPT will occur when the field brings down the triplet ipi 
and closes the gap A = g^^\0) at a critical field 



cO 



A-8 + ^A2-A2l6(A-4) 



5+(l-A2) 



(4) 



where Ag = g-/g^. Further increase of the field will 
exhaust all the singlet and so fully polarize the GS at 



Hf = 



A + 8+./A2 + A2l6(A + 4) 



fl+(l-A2) 



(5) 



which is valid for all values of A. For a weak OILSC 
< A < Ac+, the field brings out the components 
<P3,<P2, 'Po in turn at three different critical points. 
For A < 0, the single-ion energy Eq of the singlet po 

rises when the field becomes stronger. The triplet com- 
ponent p3 has lower energy in zero field but will be raised 
by the field in a quicker way than (po- The energy com- 
petition of </?3 and (po results in new kinds of QPT's. In 
the four-component GS a weak negative OILSC does not 
expel the singlet far away from the triplet. In an increas- 
ing field (^3 rises in energy more quickly and soon goes 
beyond pQ, so 953 gets out of the GS before po. There are 
three critical fields consecutively exhausting P3,PQ and 
P2 from the GS. However, a stronger negative OILSC 
leaves a longer distance between the triplet and singlet. 
Before p3 covers this distance to pa, the singlet po has 
already risen completely out of the GS, which happens 
in the first critical point. But p, soon draws very near 
to Po such that po becomes close enough to the GS, in 
which p, still lies, and is drawn back into the GS. This 
brings about the second QPT. Afterwards, the field over- 
whelms the influence of OILSC and pump out p:^, pQ and 
P2 one by one, which gives other three QPT's. Therefore 
the system undergoes five consecutive QPT's (5-QPT) 
in this case, the variations of the components in the GS 
are: 1230^123^1230 (and 1203)^120^12^1, where 
we denote the component pi by number i and each arrow 
indicates a QPT. Such unusual 5-QPT occur near Ac- 
(Asc < A < Ac-, Asc = -0.712, for g, = 2.0, gt = 1.0). 
When the OILSC is negatively stronger than Jc-, the 
singlet is gapful and does not exist in the GS before the 
field is applied. But the quicker rising of (^3 will get 
the singlet involved in the GS, as what happens in the 
second QPT of the 5-QPT. In such case, four consecu- 
tive QPT's (4-QPT) will be observed: 123^1230 (and 
1203)^120^12^1. A negatively stronger OLSC than 
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Xq will expel the singlet too far away from the GS com- 
posed of the triplet. The triplet component tps itself will 
get out of the GS before it becomes adhesive and draw the 
singlet into the GS. We then have only the two QPT's. 

We present a detailed phase diagram in Fig.l, the in- 
set illustrates the 5- and 4-QPT. 5-QPT exist if C_Q 
tends to rise up at C_, with requirement on the g fac- 
tors 5t < (1 - ^gb)gs/ (1 + Ag5) = 0.5165s, where Ags = 
3[(1 + w)ao]^/^/(47r) and w = 2/3, satisfied for most typ- 
ical spin-orbital systems. When the triplet 993 attracts 
the singlet (po into the GS, some density of ^3 themselves 
are attracted upwards a bit by <po, which causes a quicker 
increase in magnetization. This attraction-and-counter- 
attraction is more sensitive near the point Q where tps 
has small proportion in the GS. A 4-QPT magnetization 
near point Q is plotted in Fig.2B. The magnetization in 
the 5-QPT case has less apparent changes in the first two 
transition points due to heavier density of v's- 

Besides the exact singlet phase critical point (4) and 
full-polarized point (5), other critical points of the QPT 
can be obtained from Wiener-Hopf method [21] 

fjC^NQ ^ H++Tk+G\,^_,{X,H+), 

H^+^ = Ho + Tko[wGl _i 3{X,Ho) + G^_i 1 .iX,Ho)], 

2'2-2 2!2'2 

H^Q =Ho+ TkoiwGl, 1 oi\Ho) + Gl, , 3 (A, Hq)], 
where ko = 3g+^, r = 1/(477^), Ho = (ao - X/2)g+^, 
A± - 3±y(A^-3i±iA2)A2 + A2 

= 3(1T1)/25^(1_9±1A2) ' 

Gi,mA^^ H) = ^A/2 + mg+H/2 + n^X?VgilP/2, 

_ 4 3(^^3)/.^.^^^^ 
^ 3(iTi)/2g+^'^ A±-3(iTi)/2g+F±^ ' 

and A± = 8 In 2 ± A. H± and Hq are the contribu- 
tions of infinity Fermi points (FP), the other terms come 
from the revision of the finite FP. For gs = 2gi. H+ 
can be simplified as H+ = [X\ — 9A2]/[25f+A+]. The 
above analytic results are compared with the numeri- 
cal ones in Fig.l, which shows an excellent accuracy. 
For C+NQD and C+MQ, deviations only occur near 
C+ due to small FQ invalid for Wiener-Hopf method, 
but an expansion based on small FP can be carried out 
and H^ = (2c/7r)5-i(Ac+ - A)3/2, c = 1 for C+N and 
c = (ln2)/3 for C+M. Setting A = in H^^ , 
and Hf will recover the very accurate expressions of the 
SU(4) model obtained in Ref. [10]. 

A remarkable progress in experiments of spin-orbital 
systems is the observation of orbital ordering in some 
transition- metal compounds [2] , which can be also found 
in the present model. The expectations of the spin 



s and orbital angular momentum I, the magnetization 
= gsS^ + gj^ can be obtained for a singlet site 

(s-)o = -(Oo = g-H[2Xj{H/X}]-\ (Af-)o = g-{^o- 

Many transition-metal elements involve a strong OILSC 
and the GS of the system can be located in the singlet 
phase, in which all sites arc occupied by the singlet. Con- 
sequently, the spins will be aligned ferromagnetically in 
the z direction while all the orbitals are ferromagnetically 
ordered in the opposite direction. Therefore, both orbital 
ordering and spin ordering can be observed. In addition, 
unlike the conventional singlet, the magnetization {M^)o 
is no longer zero in the field, as shown by the dash lines 
of curves e and / in Fig.2A. {M^)q increases almost in a 
linear way due to a small value of H/ A before the singlet 
phase is transited to another phase at Hco- But the in- 
creasing rate d{M^)o/dH becomes smaller when the field 
gets stronger. This will lead to unsymmetric increasing of 
the magnetization in weak and in strong fields. The dot- 
ted regions of curves c, d and e in Fig.2A show a quicker 
climbing than the regions thereafter, which happens in a 
same phase with only two components (pQ and (pi in the 
GS. Curve d only involves these two components from the 
beginning and provides a full view of this unsymmetry. 
Besides the slowdown of d{M^)o/dH , another contribu- 
tion to the unsymmetry comes from the variation in the 
energy-difference increasing rate. The energy difference 
\Ei — Eo\ between the GS components ipi and (pQ varies 
in a slower rate d\Ei — Eo \ /dH when the field reaches 
higher. 

We also present two typical magnetization for nega- 
tive OILSC, shown by curves a and b in Fig.2A. The 
corresponding QPT's are: (a) 123-^12-^1; (6) 1230 (and 
1203)— > 120— >12-->1. One can notice the crossover of the 
two curves in phase 12 before full polarization, where only 
two components tpi and (p2 are left in the GS. Curve b lies 
in four-component phase 1230 in the absence of the field, 
while a stronger OILSC expels the singlet off the GS and 
curve a starts from triplet phase 123. The negative A 
makes po rising in energy level while p2 sinks when the 
field is applied. But the lowest state pi sinks fast and so 
all the other components are moving out of the GS. Thus 
there is one more component po getting out of the GS 
for case b than case a. Therefore, the proportion of the 
most energetically favorable component (pi, with a posi- 
tive magnetization g+/2, increases more quickly in case 
b than in case o, which gives a quicker rise of magneti- 
zation in curve b before the first QPT. The component 
(p3 in case b has smaller proportion in the beginning of 
a four-component phase 1230, while it has larger pro- 
portion in case a starting from a three-component phase 
123. Consequently, (pz gets out of the GS earlier in case 
6 and the first QPT occurs ahead of case a. But the rise 
of curve b soon slows down after its earlier first QPT, 
since the (po has weaker negative magnetization than (ps 
in case a. After p>3 gets out in case a, curve b soon rises 
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beyond curve a again as it has two components (p2 and 
(/?o moving out while only one component in case a. 
As a result, case h has higher magnetization than case a 
when they both come to phase 12. But in this phase, ip2 
rises away from the lower state ipi always more quickly 
in case a with a stronger A, as shown hy d\E2 — Ei \ / dH 
in Fig.2B. Therefore, 1^2 gets out faster in case a, and its 
lower magnetization rises in advance of case h and reaches 
the saturation point first. This results in the crossover. 

Finally we stress only in the context of the different g 
factors of spin and orbit that the OILSC leads to these 
properties. If g^ = gt [17], the novel phenomena including 
the five consecutive QPT's, the spin and orbtal ordering, 
nonzero magnetization of the singlet, the unsymmetric 
magnetization and the magnetization crossover, will all 
disappear. 
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External field H 

FIG. 1. Phase diagram of quantum phase transitions (QPT's) in magnetic field, = 2.0, gt = 1.0. The number i labels the 
state </5i, e.g., the state components in phase 0123 are v'oViVa'^ss in which the singlet 950 is energetically the most favorable 
whereas the triplet component 993 is the least favorable. The discrepancy between the analytic and numerical critical fields 
is not visible for most region of C+NQD and C+MQ. The boundary of singlet phase and the full-polarized point are exact. 
INSET: the five consecutive QPT's (dark dots) and four QPT's (dark boxes), their full-polarized points are relatively far away. 
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FIG. 2. (A) Typical magnetization behaviors in different QPT's, Qs = 2.0, gt ~ 1.0. The dashed hnes are gapped singlet 
phase with a nonzero magnetization , both spin and orbital ordering can be observed in this phase. Curve d provides a 
full view of unsymmetric growth of M^, the dotted- line region (phase 01) in weak field climbs faster than the solid-line region 
(phase 10) in strong field. Also note the crossover between curve a and b before full polarization. (B) Magnetization of four 
consecutive QPT's, A = —0.76, Hc4 ~ 3.65, the numbers denote the phases as in Fig.l. (C) Comparison on increasing speed of 
energy difference E2 — Ei in the field for A = —1.5 and A = —0.5, the former always has a larger speed. 



